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We investigate the temporal and spectral correlations between flux and 
anisotropy fluctuations of TeV-band cosmic rays in the light of recent data taken 
with IceCube. We find that for a conventional distribution of cosmic-ray sources 
the dipole anisotropy is higher than observed, even if source discreteness is taken 
into account. Moreover, even for a shallow distribution of galactic cosmic-ray 
sources and a reacceleration model, fluctuations arising from source discreteness 
provide a probability only of the order of 10% that the cosmic-ray anisotropy 
limits of the recent IceCube analysis are met. This probability estimate is nearly 
independent of the exact choice of source rate, but generous for a large halo size. 
The location of the intensity maximum far from the Galactic Center is naturally 
reproduced. 



Subject headings: cosmic rays 



1. Introduction 

The anisotropy in the arrival directions of cosmic ra ys in the TeV band has received re- 



newe d attention in recent y e ars. A variety of experiments (|A 
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2012[ ) have reported results on the amplitude of the first harmonic in the 



siderial anisotropy, which is typically found to be slightly less than 0.1%. All experiments 
see only part of the sky, and it is not easy to reconstruct an allsky dipole anisotropy, largely 
because the energy dependence of the acceptance is different among the various detectors. In 
any case, the true dipole anisotropy can be slightly higher than the sidereal first harmonic, 
i.e. < 0.1%. In addition, significant small-scale anisotropy was observed that has not met an 
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accepted explanation to date, although it is tem pting to attribute it to the local structure of 



the turbulent magnetic field in the Galaxy (e.g. iGiacinti fe Siglll2011t iGiacinti et al.l 120121 ). 



Here we investigate whether or not the dipole anisotropy can be reproduced with mod- 
els of galactic cosmic-ray propagation that are tuned to fit available data on secondary- 
to-primary ratios and the survival fraction of unstable isotopes. For that purpose we 
use a Green-fu n ction m ethod that solves a time-dependent diffusion equation, similar to 
Blasi &: Amatol (J2012bJ). Diffusion equations of this type are derived by averaging the orig- 
inal Fokker-Planck transport equation for the isotropic part of the cosmic-ray distribution 
function. By construction, anisotropy in this treatment arises only from the diffusive (and 
possibly drift or convective) flux and is always dipolar in nature. Modeling the small-scale 
anisotropy, on the other hand, would be described with a transport equation derived from 
the original Fokker-Planck equation for the anisotropic part of the cosmic-ray distribution 
and is not considered here. 

Our time- dependent treatment permits accounting for a discrete nature of cosmic-ray 
source. A similar technique was used to mi 



ray electrons (1 
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ray source distributions in galactocentric radius, using the appropriate transport parameters, 
and compute the expected temporal evolution of the anisotropy. This approach permits a 
realistic assessment of the likelihood to find a certain anisotropy for a specific choice of pa- 
rameters. The occasionally discussed mean variance of the anisotropy is a considerably less 
useful quantity, because the distribution of anisotropy amplitudes is highly skewed. 

We are particularly interested in the temporal and spectral correlations between flux 
and anisotropy fluctuations. Our approach is ideally suited to recover these correlations and 
investigate their effect, when we compute flux and anisotropy over a finite spectral acceptance 
and compare with observational results. 



2. Calculation of particle spectra 

2.1. Cosmic-ray ions 

The differential density of TeV-band cosmic rays obeys a continuity equation that can 
be written as ^ ^ ^ 

dt r 2 dr \ dr j ^ ^ ^ Anr 2 ^ ^ 
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Here, D = cA m f p /3 is the spatial diffusion coefficient, r is the distance from the source, and 
Q(t, E) is the differential production rate at the source. Note that we neglect stochastic 
reacceleration, advection, and energy losses. The first and the last of that list should only be 
important at GeV energies and below, at least for cosmic-ray protons. The case of advection 
is less clear, but probably also less important than it may be in the GeV band on account 
of the energy dependence of diffusive propagation. 



Using standard methods (jKardashevlll962l ; iGinzburg fc Syrovatskiilll964j ) one derives a 
general solution to the transport equation, assuming impulsive injection of particle at time 
t = 0, 



N(r,t,E) 



e(t) 



Q{E) exp 



(2) 



If more than one source contributes to the local cosmic-ray flux at any time, their individual 
contributions must be calculated using equation [2] and then summed. 

We now modify the solution to account for a finite size and lifetime of sources and for 
escape from the Galaxy. It is known that assumi ng a point sou rce in space and time can 



lead to singularities in the cosmic-ray density (e.g. lMertschll201l[ ). Real cosmic-ray sources, 
whatever their nature, have a finite spatial extent and confine the accelerated particles for 
some time, for the acceleration is not instantaneous. The most simple method to represent 
extended sources is to assume that they have the shape of a thin, spherical shell of radius 
R. The correspondence to the shell of a young supernova remnant is obvious, but the 
applicability of the ansatz is not limited to SNR. In any case, averaging over the shell yields 



N(r,t,E) 



8(t) 



(4tt D t) 
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Q(E) 
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(3) 



Note that this solution still has a singularity for t — > at r = R, which will disappear when 
we average over a finite lifetime, because asymptotically N oc t~ 05 for £—>•(). 

We do need to account for source size in computing the anisotropy, which is given by 
the diffusive flux, 
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coth 



rR 
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(4) 



For small argument of the coth term, the first and last term on the RHS cancel each other. 

For the integration over the time of source activity, we choose the most simple form 
of a temporal profile, namely a constant source rate over an age interval [Ti;T2] with the 
constraint that both T% and T2 be positive. The time integral can be easily performed 
analytically, but leads to very small differences of exponentials and error functions that are 
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more difficult to numerically implement than a simple numerical integration. A fact that 
we shall employ below is that either escape is irrelevant or the effects of the finite extent of 
lifetime of the sources are negligible, provided the source lifetime is much smaller than the 
escape time. 

Diffusive escape from the Galaxy can be modeled using Dirichlet boundary conditions at 
some distance from the midplane of the Galaxy, N(z = ±H) = 0. Obviously, the boundary 
conditions break the spherical symmetry of the p roblem. Whereas for very rare sources 
analytical approximations may ease the treatment (jPohl fc Eichlerl 120111 ). in the TeV band 
we expect cosmic-ray sources to appear with a rate ~ 1/(100 yr), implying a relatively small 
distance to the nearest source. Consequently, we attempt a full solution by using the method 
of mirror sources. Since we have two boundaries, an infinite series of mirror sources with 



alternating sign is needed (IBlasi fc Amatol l2012al ). 



For the mirror sources only the Gaussian factor in equation [3] is relevant. It is accurate 
to second order in the ratio of source lifetime and escape time to use an average age, r = 
(T 2 + Ti)/2. Mirror sources modify the z-dependent Gaussian factor in Equation [3] as 
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which, as is described in Appendix |Aj can be approximated with good accuracy as 
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2.2. Cosmic-ray electrons 



It is well known that the flux and spe ctrum of cosmic-ray electrons from discre te sources 
are particularly sensitive to fluctuations (jPohl fc Espositolll998l ; iPohl et al.ll2003l ). but the 
consequences for the electrons anisotropy are rarely discussed in the literature. In fact, the 
anisotropy in the intensity of cos mic-ray electrons at ener gies between 60 GeV and 480 GeV 
has only recently been measured (lAckermann et al.ll2010l ). In the case of electrons, we need 
to account for their energy losses, but may neglect escape, at least above ~60 GeV. Hence, we 
may describe the propagation of electrons at energies higher than a few GeV by a simplified, 
time-dependent transport equation 

_£ _ ^L {hE 2 N e ) — D E a V 2 N e = Q e (7) 

where we consider continous energy losses by synchrotron radiation and inverse Compton 
scattering, a diffusion coefficient D E a dependent on energy, and a source term, Q e . Green's 
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function for this problem can be found in the literature (IGinzburg fc Syrovatskiil 119641 ) . 
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As in the case of cosmic-ray ions we assume the individual sources inject a spectrum 
oc E~ s at constant rate for a defined period of time, and so the average over age from T to 
T 2 is 



N e = q e E~ 
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and r is the distance between source and observer. 

As in the case of ions, we can easily average over a shell of radius R and obtain 
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A general analytical solution is very difficult to obtain. However, an explicit integration is 
unnecessary if (T 2 — T\) <C T, provided the energy-dependent upper limit of integration is 
observed. If s > 2, as is the case for cosmic-ray electrons, the integrand approaches zero at 
t (bE)^ 1 . Then it is accurate to second order in (T 2 — T\)/Ti to set 
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where r = (T 2 + Ti)/2. The energy-loss time of electrons at 500 GeV in a 10-/iG magnetic 
field is about 3 • 10 5 years, much longer than the source life time T 2 — T\ that we consider 
here. If bEt is a small parameter, we can expand the integrand in Equation [12j Noting that 
for bEt < 1 

A~D E a t, (14) 

we find that Equation [12] assumes the same form as the time integral over Equation GU the 
corresponding expression for cosmic-ray ions. 
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3. Intermittency and anisotropy 



We now run 100 simulations for each parameter combination. Unless the number of 
cosmic-ray sources is very small, and hence intermittency effects dominate, the source dis- 
tribution is a critical factor in anisotropy studies. Often, a function of the type 



Q{ r Gc) 



R<7) 



exp 
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R r . 



(15) 



is used to represent the density distribution in galactoce ntric radius, rp c- The parameters a 
and b derived by fitting obser ved distributions of pulsars (jLorimerll2004j ). supernova remnants 
( jCase fc Bhattacharyalll998l ). or other perceived tracers of cosmic-ray acceleration. 



It has been noted earlier that the observed intensity of diffuse Galactic gamma rays 
appears to require a substantially flatter source distribution ( IStrong et al.ll2000l ). although a 
metallicity gradient, and hence a gradient in the Xgo conve rsion factor of CO line intensity, 
could explain the ga mma-ray data as well (IStrong et al.ll2004j) . The comprehensive analysis of 
Fermi-LAT data bvlAckermann et al. also suggests a gradient in Xc p, albeit a weaker 



trend as in IStrong et al.l (120041 ) . A similar conclusion was reached by lAckermann et al. 



Here we select from the literature three models of the cosmic-ray source distribution 
i n the Galaxy. Our sample comprises the relatively flat dist r ibutio n used by IStrong et al. 



( 120001 ). the S NR dis t ributi on derived by lCase fc Bhattacharya! ( 119981 ). w hich is slightly flatte r 
than that of iGreenl (I2012[). and the puls ar distribution as presented by iTrotta et al.l ( 1201 ll ). 
which is close to that of lLorimerl ( 120041 ) . Figure [T] permits a visual comparison of the three 
models. 



According to ITrotta et al.l ( 120111 ). which we choose for the definition of standard pa- 
rameters, we have well-constrained propagation parameters, such as the injection index 
s = 2.4 ± 0.1, and for the diffusion coefficient 



D = D 



4 GeV/ 



(16) 



with 5 = 0.3 ±0.05 in reacceleration models, whereas a degeneracy arises with the halo size. 
In units of kpc, 3 < H < 9, and 



D ~ (1.2 + 1.3 if) • 10 28 cm 2 s 



(17) 



within about 20% error margin. We try 4 values of H and derive the 95%-confidence minimum 
of anisotropy in the 10-100 TeV band. Standard parameters are summarized in TableUJ Note 
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Fig. 1. — Distribution of the cosmic-ray source density in galactocentric radius, according 
to 3 estimates used in the literature. All curves are normalized to unity at the solar circle. 



Table 1. Standard parameters, mostly taken from the study of [Trotta et al.l (1201 if ) 



Parameters 


Symbol 


Value 


Injection index 


s 


2.4 


Energy dependence of diffusion 


6 


0.3 


Source distribution 


a 


1.25 


Source distribution 


b 


3.56 


Source rate 


Pq 


10~ 2 yr" 1 


Source lifetime 




2 ■ 10 3 yr 


Source radius 


R 


10 pc 
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that these parameters were determined using the sharply peaked source distribution labeled 
Trotta et al. in Figured! The anisotropy and its variation is shown as function of the particle 
energy in Figure [U full results are detailed in Table [2J To be noted from Figure [2] is the 
skewness in the distribution of anisotropy levels, which renders the use of a mean fluctuation 
amplitude in the anisotropy less than useful. 



It is obvious that using the viable propagation parameters of iTrotta et al.l (120111 ). the 
limits on anisotropy may be met only if either the TeV-band cosmic rays are relatively heavy 
or the source distribution is very flat. A shallow source distribution is certainly possible, 
but was not the one used to find the propagation parameters best suited to reproduce B/C 
or 10 Be data. 



3.1. Shallow source distribution 



To verify the reliability of our findings, we test with results of iPutze et al.l (120101 ) who 
used a different propagation model and a constant source density. We do not use their 
reacceleration/convection model, because in 2-D geometry with radial boundary the best-fit 
parameters indicate that advection has little importance above a TeV. 

For the reacceleration model the scaling of the diffusion coefficient with halo size (in 
kpc) is 



D ~ H 1.04 • 10 



28 cmV 1 



with 1. 5 < H < 6. Num erically, the scale factor of the diffusion coefficient is similar to 



that of ITrotta et al.l ( 1201 ll ). but at 20 TeV the diffusion coefficient is actually a factor 2.6 
smaller on account of the weaker energy dependence of diffusion. O ther paramet e rs are 
summarized in Table [3j We use the relatively flat source distribution of IStrong et al.l (120001 ) 
out to tqq = 20 kpc, and henceforth refer to it as the shallow source distribution. 

Results are summarized in Table SJ To be noted from the table is the strong reduction 
in anisotropy resulting from the moderate (factor 2.6) decrease in the diffusion coefficient 
(cf. Tabled in particular the runs for Strong (2000)). 



1 10 100 1000 

Energy in TeV 



Fig. 2. — Anisotropy for standard parameters (Trotta et al.) and H = 5 kpc. The red band 
indicates the central 90% of the fluctuation range. The solid line marks the median of the 
distribution, the dotted line describes the anisotropy amplitude in a randomly selected run. 



Table 2. Proton anisotropy at 20 TeV total energy for various parameters and 4 choices of 





halo size, H 


in kpc. 








Parameters 


Confidence 


H=3 


H=5 


H=7 


H=9 


Standard 


95% lower limit 


0.29% 


0.34% 


0.55% 


0.46% 


Standard 


Median 


0.60% 


0.72% 


0.81% 


0.95% 


S = 0.25 


95% lower limit 


0.18% 


0.19% 


0.35% 


0.35% 


S = 0.25 


Median 


0.40% 


0.48% 


0.54% 


0.60% 


Case & Bhattacharya 


95% lower limit 


0.14% 


0.24% 


0.23% 


0.29% 


Case & Bhattacharya 


Median 


0.44% 


0.55% 


0.62% 


0.70% 


Strong (2000) 


95% lower limit 


0.087% 


0.07% 


0.14% 


0.12% 


Strong (2000) 


Median 


0.38% 


0.35% 


0.37% 


0.43% 
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Table 3. Parameters for a shallow source distribution, mostly taken from Putze et al. 



Parameters 


Symbol 


Value 


Injection index 


s 


2.415 


Energy dependence of diffusion 


5 


0.235 


Source rate 


Pq 


10~ 2 yr" 1 


Source lifetime 




2 ■ 10 3 yr 


Source radius 


R 


10 pc 
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Fig. 3. — Anisotropy for a shallow source distribution with parameters as listed in Table [3] 
and halo size H = 4.5 kpc. As in Figure |2} the red band indicates the central 90% of 
the fluctuation range, the solid line marks the median of the distribution, the dotted line 
describes the anisotropy amplitude in a randomly selected run. 



Table 4. Proton anisotropy at 20 TeV total energy for a shallow source distribution. 



Parameters Confidence H=1.5 H=3 H=4.5 H=6 



Standard 
Standard 
Source rate 
P Q = 1CT 3 yr- 
Source rate 
P Q = 10- 4 yr" 



95% lower limit 
Median 

95% lower limit 
Median 

95% lower limit 
Median 



0.026% 0.023% 

0.12% 0.10% 

0.022% 0.037% 

0.12% 0.13% 

0.034% 0.054% 

0.13% 0.15% 



0.032% 0.043% 

0.14% 0.14% 

0.03% 0.035% 

0.12% 0.16% 

0.035% 0.054% 

0.13% 0.16% 




Fig. 4. — Light curves of proton flux at 2 energies (top panel), proton anisotropy (middle 
panel), and electron anisotropy at two lower energies (bottom panel). The black bars and 
dotted horizontal lines indicate the median of the anisotropies. 
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4. Light curves 

4.1. General behavior 

Having established that a rather flat source distribution is necessary to explain the low 
anisotropy amplitude measured at multi-TeV energies, we now discuss light curves. We want 
to infer the temporal and spectral correlations between flux and anisotropy fluctuations, 
which must be taken into account when computing the flux and anisotropy over a finite 
spectral acceptance and comparing with observational results. 

A number of insights can be gained by first discussing light curves at specific ener- 
gies. Figure H] shows the temporal variations of proton flux, proton anisotropy, and electron 
anisotropy for the standard parameters with halo size H = 3 kpc and a shallow source 
distribution (see Table [3]). Our findings can be summarized as follows: 



• The proton flux displays s mall-amplitude fluctua tions with rare, short-lived spikes, a 



result already discussed by iBusching et al.l (120051 ) 



• Spikes in flux are generally accompagnied by spikes in anisotropy at all energies and 
reflect the appearance with nearby, young cosmic-ray sources. The events may change 
the diffusion coefficient on account of streaming instabilities, in which case this calcu- 
lation would overestimate the anisotropy during spikes. 

• Deep, short-lived dips can be observed in the anisotropy, which have no counterpart in 
the flux and do not appear at all energies. These dips arise from a random balancing 
of the contributions of the dominant cosmic-ray sources and are therefore a hallmark 
of discrete sources. 



1.8 TeV 




-5 -4 -3 -2 -1 
time in Myr 



Fig. 5. — Distribution in Galactic longitude of the direction of maximum proton intensity. 
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• There is little correlation between the small-amplitude fluctuations in the anisotropy 
at different energies, which implies that care must be exercized when comparing with 
airshower data, which carry contributions from cosmic-rays over a wide range of ener- 
gies. 



• The p redicted electron anisotropy is lower than the upper limits published by lAckermann et al 



(120 lOl ). but higher by a factor ~5 than the GALPROP predictions shown in their Fig- 
ure 9. GALPROP does not properly account for discrete sources and therefore would 
underpredict the anisotropy. 



The direction of maximum intensity is almost always in the Galactic plane, because 
the distribution of cosmic-ray sources with the vertical coordinate, z, is very narrow. In 
Figure [5] we show how the anisotropy direction varies with time, using the same simulation 
as for Figure H] for ease of comparison. To be noted from the figure is the correspondence of 
the deep dip in anisotropy at -3.2 Myr with a 360-degree swing of the anisotropy direction. 
Generally, the variations in the anisotropy direction are large and in general not achromatic, 
although the difference in directions at 1.8 TeV and 110 TeV rarely exceeds 30 degrees. 

Ground-based observatories do not see the entire sky and therefore resort to measuring 
the first harmonic in right ascension, while integrating over the declination band that is 
accessible to them. The projection of a first harmonic in galactic longitude depends on the 
direction of maximum intensity and on the declination band in question. 



4.2. Comparison with IceCube results 

The IceCube col laboration has pres ented cosmic-ray anisotropy measurements in two 



broad energy bands (lAbbasi et al.ll2012l ). based on a total of about 30 billion cosmic-ray- 
induced muon events. Airshowers caused by heavy primaries approximately behave like the 
superposition of airshowers produced by the individual nucleons. Therefore, the detection 
with IceCube of a muon with a certain energy requires a heavy primary of much higher 
energy than it would, if the primary were a proton. In terms of propagation, heavy cosmic 
rays with charge number Z and energy E will behave like protons of energy E/Z. Thus, 
light and heavy cosmic-ray nuclei that have similar propagation properties, and hence similar 
anisotropy, will also account for muons of similar energy in IceCube. We have seen in Figured 
that a factor 2 in energy on average yields a 20-% change in anisotropy. We have also seen 
that the spectrum of cosmic rays shows little variation, except possibly during intense spikes 
in flux. It is therefore a good approximation to only consider cosmic-ray protons. 
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1 10 100 1000 

E in TeV 



Fig. 6. — The distribution in energy of the primary particle for the high-energy (red line) 
and low-energy (blue line) event selections, derived under the assumption that all primary 
cosmic rays are protons. 




Fig. 7. — Top panel: Expected sidereal anisotropy in IceCube data at declination 6 C = —45° 
for the high-energy selection (in red) and low-energy band (in blue). The bars on the right 
indicate measured limits. The source rate is assumed to be one per thousand years. Bottom 
panel: Location in Right Ascension of the maximum intensity in IceCube data for the high- 
energy selection (in red) and low-energy band (in blue). 
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The IceCube collaboration was kind enough to provide us with the distribution in energy 
of the primary particle for their 2 muon event selections, derived under the assumption 
that all primary cosmic rays are protons. These distributions are also shown in Figure El 
Obviously t he median pr i mary particle energies are lower than for the mixed composition 



assumed in lAbbasi et al.l ( 120121 ). but these are the distributions we need to combine with 



our proton-propagation calculation described in the preceding sections. 

IceCube reports the sidereal first harmonic in the cosmic-ray intensity averaged over 
the declination range —25° < 9 < —72°. Obviously the dependence of acceptance on 
the primary energy is not constant over this declination range, but we will ignore that 
and project the cosmic-ray dipole anisotropy onto the sidereal band at the characteristic 
declination 6 C = —45°, which is approximately the median declination of events accepted in 
the IceCube analysis. Note that we compute the direction of the cosmic-ray dipole anisotropy 
(cf. Figure [5]), which displays considerable variability, and so does its projection onto the 
sidereal band. 

A dipole anisotropy of amplitude 5 pointing at declination 6 d provides a sideral first 
harmonic at declination 9 C of amplitude 

6 cosflrf cosflc 

Voj(0 c ) = . — 19 

1 + o sin v d sin t) c 

The denominator can essentially always be set to unity. We simulate the cosmic-ray intensity 
for 50 Myr to obtain a statistically meaningful data set. In Figure [7] we show the expected 
anisotropy in the IceCube data over 20 Myr only for a source rate Pq = 3 • 10~ 3 yr~ 4 . 

The IceCube collaboration reports the following sidereal anisotropies: 

5 obs (low energy) = (7.9 ± 0.1 stat ± 0.3 sys ) ■ 10~ 4 

5 obs (high energy) = (3.7 ± 0.7 stat ± 0.7 sys ) • 10~ 4 (20) 

The expected anisotropy is simultaneously below both values for 8.2% of the simulated time 
period of 50 Myr. Even if we account for the uncertainties and set the anisotropies to 8.2- 10~ 4 
and 4.8 • 10 -4 , respectively, implying that individually there is only a 16% probability that 
the true anisotropy is higher, they are met at only 12.9% of the time. 

Also shown in Figure [7] is the location of peak intensity in the sidereal band studied by 
IceCube. To be noted from the figure is that in particular during anisotropy lows the peak 
intensity is found far from the Galactic-Center direction at RA~ 17. 8h and not necessarily 
coincident for the 2 energy bands. Both results reproduce the IceCube observations. 

To verify that this results does not significantly depend on the assumed source rate 
and halo size, we have varied the rate between Pq = 3 ■ 10~ 4 yr _1 and Pq = 3 ■ 10 -2 yr~ 4 , 
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Table 5. Summary of anisotropy results. The first two columns give the assumed source 
rate and halo size. The next three columns give the probability, that the anisotropy is 
below that measured with IceCube, and the anisotropy median at low energy and high 
energy. The last column identifies the corresponding lines in Figure |HJ 



Rate Halo Prob. Med. LE Med. HE Line 
[K^Myr- 1 ] [kpc] [%] [1(T 4 ] [10~ 4 ] 



0.3 3 8.2 7.1 13.5 Black 

0.3 7 3.3 8.3 15.7 Blue dotted 

3 7 4.2 7.7 15.2 Red dotted 

3 3 8.2 6.8 12.8 Red 

1 3 10.4 6.8 13.2 Blue 

10 3 9.1 6.2 12.0 Green 




P(<A) 



Fig. 8. — The integral probability, P(< A), to measure an anisotropy lower than a certain 
threshold, A. Solid and dotted lines refer to the low-energy data selection of IceCube, dashed 
and dot-dashed line are for the high-energy data. The parameters of each curve are listed in 
table E 
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and the halo size between H = 3 kpc and H = 7 kpc. In no case was the probability 
to simultaneously meet both anisotropy values higher than 11%. Table [5] summarizes the 
results, and Figure E] displays the integral probability to measure an anisotropy lower than 
a certain threshold, A. To be noted from the figure is that varying the source rate does not 
change the probability distribution in a significant manner, but choosing a large halo size 
appears to reduce the probability of a low anisotropy. 

As last test, we repeated one simulation (H = 3 kpc and Pq = 3 • 10~ 3 yr _1 ) for 
a perfectly flat source distribution out to r GC = 20 kpc. The probability of meeting the 
Ic eCube upper l i mits w as with 11% only marginally er than for the shallow source distribution 
of lstrong etaD feoooh . 



5. Discussion 



This is not the first calculation of cosmic-ray anisotropy in the TeV band. The main 
advantage of this study lies in the proper accounting for the correlations between flux and 
anisotropy variations, a restriction to cosmic-ray propagation parameters that were deemed 
viable in dedicated studies, and an accurate folding with the acceptance distribution of a 
recent observational study. 

We have specifi cally compared with recent measurements with the IceCube experiments 
(lAbbasi et al.ll2012l ). which leads us to the following conclusions: 



• Consistent with earlier findings, a cosmic-ray source distribution in the Galaxy that 
fits to the deduced populations of pulsars or SNR, combined with propagation param- 
eters that reproduce the observed secondary-to-primary ratios, will lead to a dipole 
anisotropy that is higher than observed, even if source discreteness is taken into ac- 
count. 

• Also in agreement with earlier studies, a shallow energy dependence of diffusion is 
necessary, thus requiring a reacceleration model to fit the GeV-band Boron-to-Carbon 
ratio. 

• For a shallow source distribution, fluctuations arising from source discreteness provide 
a probability of the order of 10% that the cosmic-ray anisotropy limits in the high- 
energy and low-energy bands of the recent IceCube analysis can be simultaneously 
met, approximately independent of the exact choice of source rate, but less likely so 
for a large halo size. 
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• The location of the intensity maximum far from the Galactic Center and at different 
locations for the 2 energy bands are naturally reproduced when source discreteness is 
taken into account. 



It is unclear what type of cosmic-ray source in the Galaxy would have a distribution 
in the Galaxy as shallow as required in our study. The question arises whether or not 
modifications to the description of diffusive transport can render a steep source distribution 
viable. The critical constraint appears to be imposed by very low dipole anisotropy in the 
high-energy band of IceCube data, i.e. at approximately 100 TeV energy per nucleon. Is 
it possible that the extrapolation to this energy of the cosmic-ray propagation parameters 
does not hold? 



Recent data suggest th at cosmic-ray spectra (Ahn et al.ll2010i ; lAdriani et al.ll201ll ). and 
the Boron-to-Carbon ratio ( lObermeier et all 120121 ) . flatten above a few hundred GeV/nuc. 
At this time we do not have a prevalent interpretation of both the cosmic-ray spectral hard- 
ening and the flat B/C ratio above a few hundred GeV/nuc. The spectral hardening of 
cosmic-ray primaries could simply be a fluctuation arising from source discreteness, but in 
that case the Boron-to-Carbon ratio s hould fall off more stee ply, because the flux of cosmic- 
ray secondaries fluctuates very little (IBusching et all 120051 ) . S pectral hardening can also 



reflect variations in t he production sp ectra of different sources (IBrecher fc Burbidgd 11972 



Biisching et ajj |2Q0l| ; lYuan et all 120111 ). in which case the Boron-to-Carbon ratio would be 
unaffected. If the upturn in the Boron-to-Carbon ratio is, on the other hand, due to the 
fragmentation of primaries and the re-acceleration of secondaries inside cosmic-ray sources 
(IMertsch fc Sarkarl 120091 ; iTomassetti fc Donatol 120121 ). then the energy dependence of the 
diffusion coefficient cannot be argued to be very shallow. The only possibility that would 
imply a very low anisotropy is a change of the energy dependence of t he diffusion coef- 
ficien t at a few hundred GeV/n uc, e.g. from 5 = 0.25 to 5 = 0.15 (IVladimirov et al. 



20121 ; ITomassetti fc Donatoll2012l ). Speculative though that may be, the extrapolation using 
5 = 0.15 to 10 17 eV may render viable t he notion that a signifi cant fraction of the sub-ankle 



cosmic-ray flux is galactic in origin (e.g. iPohl fc Eichlerl 1201 lh 



Recently, lEvoli et all (120121 ) suggested a common solution to the anisotropy and cosmic- 
ray gradient problems, that involves a strong toroidal magnetic field in the Galaxy. Diffusion 
perpendicular to the azimuthal magnetic field (with coefficient D± oc E 0M , as opposed to 
D11 oc E 3 ) would be the dominant transport mechanism, assuming field-line wandering and 
large-scale velocity turbulence can be ignored. The authors further posit that -D_l(?"gc) 
l/D||(rGc) an d D± oc Q(r GC ) T , where Q(r GC ) is the density of cosmic-ray sources and 
r ~ 0.85. At least up to 1 TeV particle energy, the predicted radial anisotropy appears to 
be commensurate with measured values. 
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We posit that the scenario would do little, if anything, to resolve the 3-D anisotropy 
problem for discrete sources (the authors only discuss anisotropy for continuous sources in 
a 2-D scenario). In the direction of fast diffusion,i.e. along the large-scale toroidal magnetic 
field, we see particles of the same age from sources at larger distance, thus increasing their 
contribution to the anisotropy. The volume sampled, and thus the number of sources con- 
tributing to the local cosmic-ray flux, will also increase, thus reducing the likelihood of a 
fluctuation toward low anisotropy, but not ef ficiently enough to overcome anisotropy arising 
from source discreteness. In the scenario of lEvoli et al.l ( 120121 ) we therefore expect signifi- 
cantly more anisotropy in direction of the ordered magnetic field than in radial direction. 
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acceptance distribution for their anisotropy analysis. MP acknowledges support by the 
Helmholtz Alliance for Astroparticle Physics, HAP, funded by the Initiative and Networking 
Fund of the Helmholtz Association. DE acknowledges support from the Israel-U.S. Binational 
Science Foundation, the Israeli Science Foundation, and the Joan and Robert Arnow Chair 
of Theoretical Astrophysics. 



A. Analytical appoximation of escape 



Mirror sources modify the z-dependent term as 
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that can be rewritten as 
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The solution is inconvenient, because one needs to sum over many terms to find convergence 
for all t. We are interested in the intermittency behaviour and thus need to account for many 
million cosmic-ray sources in the Galaxy, which makes computational efficiency imperative. 
We can expand the exponential in the series to give 
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where 



2nH - z s (1 



2Dt 



(A4) 



The last factor in brackets is essentially unity, unless H 2 ^> D r in which case the exponential 
yields zero, because its argument is always at least a factor H/z larger as the z-dependent 
terms in the last factor. 



We may now expand the remaining exponential for small z s to find 
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where summation over ±n removes the first-order term and otherwise yields a factor 2. 
Hence 



G ~ exp 
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As for z, the second-order term in z s can become relevant only when the exponential is 
practically zero and the term in brackets unity. Therefore, it can be neglected. A nice 
analytical approximation for the term in brackets is given by 



G ~ exp 



{Z - Z a )' 

ADt 



[1 + 1.5 a;) 1 - 25 exp [-(1.5 a;) - 97 ] 



x = 



2Dt 



(AT) 



which deviates from the true solution by less than 7.2% for x < 8, for which the term in 
brackets is down to 3.6 • 10~ 4 . 
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